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Introduction 


We  had  anticipated  that  the  39  month  period  July  1979  -  September 
1982  would  be  primarily  devoted  to  the  analysis  of  the  decay  of  energetic 
particle  events,  using  first  the  JHU/APL  IMP  7/8  CPME  (Charged  Particle 
Measurement  Experiment)  measurements,  and  then  combining  them  with  LECP  (Low 
Energy  Charged  Particle)  measurements  on  the  VGR  1/2  spacecraft,  launched  in 
August  1977.  However,  continuing  research  by  our  own  and  other  groups  reveal¬ 
ed  new  complexities  in  the  behavior  of  low  energy  «  1  MeV)  protons,  the  only 
protons  for  which  we  had  anisotropy  measurements  from  the  CPME.  Moreover,  a 
full  set  of  VGR  particle,  plasma  and  field  data  was  not  made  available  to  us 
until  last  year. 

Consequently,  we  modified  our  research  program,  concentrating  during 
the  early  period  on  subjects  relevant  to  the  problem  of  decay  phase  predic¬ 
tion.  This  background  proved  beneficial  when  we  were  able  to  work  on  the  main 
problem  this  past  year.  Some  15  papers  were  published,  and  7  others  are  cur¬ 
rently  pending  publication.  They  are  listed  (by  number)  in  the  following 
section. 


Briefly,  our  publications  fall  into  four  general  categories: 

I.  Interplanetary  Particle  Propagation:  Papers  1,  2,  8,  12,  16 

and  1 7 . 

II.  Solar  Wind  and  Geomagnetic  Storms:  Papers  3,  4,  6,  7,  9  and 
13. 

III.  Solar  Flares  and  Activity:  Papers  5,  10,  14,  15  and  21. 

IV.  Diagnostics  of  Particle  Propagation  Using  Energetic  Ions  Up¬ 
stream  of  the  Earth's  Bow  Shock:  Papers  11,  18,  19,  20  and  22. 

In  the  remainder  of  this  Report,  we  summarize  the  research  to  date 
on  the  main  problem,  that  of  predicting  the  decay  phase  of  solar  proton  events 
at  energies  >  1  MeV. 
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Theoretical  Background 

What  Information  is  contained  in  the  particle  flux  anisotropy? 
Which  component  of  the  anisotropy  gives  the  clearest  indication  of  the 
evolution  of  the  particle  event?  We  turn  to  a  global  formulation  of  the 
propagation  equations  for  the  answers,  because  the  evolution  of  an  event  is 
simply  the  competition  among  injection,  escape,  acceleration  and  deceler¬ 
ation.  The  equation  which  inter-re lates  all  these  processes  is  the  equation 
of  continuity  (conservation  of  the  particle  density  V  at  kinetic  energy  T  in 
the  presence  of  streaming  _S) .  If  we  write  this  equation  down  in  a  local  frame 
in  which  the  electric  field  is  zero,  it  takes  the  familiar  form 


it v' +  •  r  -  o 

One  such  local  frame  is  that  moving  with  the  electric  field  drift  velocity 
o 

(c/B  )  E_  x  _B .  Here  _E_  is  in  the  rest  (spacecraft)  frame,  so  E  *  -(v_  x  jl)/c  and 
consequently  the  drift  velocity  is  V_^  ,  the  component  of  the  solar  wind 
perpendicular  to  the  Interplanetary  magnetic  field.  We  have  established,  in 
Papers  (2)  and  (8),  that  for  protons  with  T  <  1  MeV  the  streaming  in  the 
spacecraft  frame  is  given  very  accurately  by 


(2) 


where  C  is  the  Compton-Getting  factor  for  non-relativistic  particles 
2(  y+  l)/3,  (~y)  being  the  slope  of  the  spectrum  differential  in  energy. 
Omitted  from  (2)  is  the  gradient  anisotropy  term,  proportional  to 
_B  x  VU,  since  it  scales  as  the  gyroradlus  and  hence  it  is  usually  unimportant 
at  non-relat  Ivist  ic  energies  (as  was  borne  out  by  our  measurments) .  By  the 
linearized  transformation  for  weak  anisotropies,  the  streaming  S'  which 
appears  in  (1)  is  related  to  in  (2)  by 
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(3) 


S'  -  S.  -  C  V  u 


which  implies  JS'  *  j>( .  We  leave  _S_  unspecified;  it  could  be  caused  by 
diffusion,  or  it  could  be  dominated  by  scatter-free  propagation.  Equation  (1) 
is  valid  in  either  extreme. 

One  way  to  see  the  physical  information  contained  in  S  is  to 

H 

integrate  (1)  over  the  volume  of  a  magnetic  flux  tube  whose  ends  are  formed  by 
its  intersections  with  the  spherical  surfaces  r  =  rQ  (the  upper  corona,  for 
example)  and  v  (the  spacecraft).  Let  these  end  areas  of  the  tube  be  Aq  and  A, 
respectively.  Then,  applying  the  divergence  theorem  to  V'  •  S'  and  recall¬ 
ing  S'  is  parallel  to  the  sides  of  the  flux  tube,  we  obtain 
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(4) 


This  says  that  he  radial  component  of  the  parallel  streaming  S|r 

(which  is  not  the  same  as  the  parallel  component  of  the  radial  streaming  S  ), 

^  II 

responds  to  the  (radial)  injection  streaming  SQ  (which,  being  the  net  flow 
across  the  surface  r  -  r  ,  contains  not  only  the  source  strength  of  new 
particles,  but  the  flow  of  previously  injected  particles  from  interplanetary 
space  back  into  the  corona),  and  the  co-moving  time  rate  of  change 
3U'/3t,  integrated  over  the  volume  of  the  flux  tube.  This  latter  integral 
is  not  3N/3t,  the  time  derivative  of  the  total  number  of  particles  in  the  flux 
tube,  because  the  derivative  is  taken  in  a  different  local  frame  at  each  point 
in  space.  The  relationship  between  3U'/3t  and  3U/3t,  the  time-rate-of-change 
of  particle  density  observed  at  the  spacecraft,  can  be  derived  from  the  equa¬ 
tions  Roelof  published  in  the  Proceedings  of  the  14th  International  Cosmic  Ray 
Conference,  (Munich,  1716,  1975),  under  the  assumption  of  negligible 

transverse  diffusion,  as  validated  in  Paper  (2): 


3U' 
3t  ' 


1U.+ 

3t 


VU  +  (V  •  V  )  CU 


(5) 
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The  three  terms  should  be  familiar.  The  first,  integrated  over  any 
volume,  does  give  3N/3t  in  the  rest  frame.  The  second  is  the  apparent  temp¬ 
oral  change  due  to  _K  x  B  convection  of  the  particle  population  in  the  presence 
of  a  density  gradient  in  radius,  longitude,  or  even  latitude,  since _V  always 

has  a  latitudinal  component  whenever  ^  does.  For  an  ideal  Archimedean  spiral 

o 

field,  V  «  r  and  is  nearly  azimuthal  well  inside  1  AU,  while  V  *  V  and  is 

2  ^ 

nearly  radial  well  outside  1  AU.  Likewise,  7  •  V_  =  4  ft  r/V  in  the  inner 
heliosphere  and  7  •  V.  *  2  V/r  in  the  outer  heliosphere,  reaching  a  maximum 
value  of  about  3  ft/2  near  1  AU,  where  ftr  *  V;  ft  is  the  sun's  sidereal  rotation 
rate.  If  there  is  interplanetary  acceleration,  for  instance  in  the  vicinity 
of  the  flare-associated  shock  wave,  then  a  fourth  (and  positive)  term  should 
be  added  to  (5). 


We  can  now  see  what  goes  into  the  parallel  velocity  anisotropy 
3  S  /Uv.  Substituting  (5)  into  (4),  we  find  that 


'o  -»*-  -  .-OX  /  *  *  I  -ft*  h  ■  m  *  4  ’  •  V  01,1  <« 


For  a  flux  tube  which  is  not  too  deformed  from  a  spiral, 
Aq/A  *  rQz/r  ,  so  the  source  Injection  term  is  important  at  1  AU  only  when 
there  is  an  initial  negative  radial  gradient,  as  strong  as  U  r~2  for  in_ 

stance.  Later  on,  as  the  gradients  weaken,  £  Is  determined  by  the  volume 
Integral,  l.e.,  (-  3U'/3t)  acts  as  an  effective  "source"  for  the  parallel 

streaming  S  r  at  radius  r.  Conversely,  if  we  took  rQ  *  1  AU  and  r  >  1 
AU,  3U'/3t,  if  positive  over  the  outer  flux  tube,  would  drive  the  streaming  at 
1  AU  inward  ( CQ  <  0),  since  the  radial  component  of  £;  becomes  small  at  large 
r  (because  Bf/B  «  I). 


Use  of  Anisotropy  Measurements  In  Predicting  the  Evolution  of  So lar  Proton 

Events 

The  full  complexity  of  Equation  (6)  does  not  always  have  to  be 

considered  In  deciding  the  nature  of  a  particle  event.  For  Instance,  the 

rising  phase  of  a  new  injection  (even  in  the  presence  of  a  pre-existing 

population)  is  easily  identifiable  by  the  strong  anisotropies  which  accompany 
the  rise-to-maximum;  no  knowledge  of  global  gradients  is  required  for  this 
signature.  Anisotropy  measurements  can  also  be  used  by  themselves  to  dis¬ 
tinguish  slowly  rising  solar  events  from  "corotating"  populations  from  a 
previous  injection.  For  example,  if  the  decaying  population  from  an  eastern 
flare  (whose  maximum  spatial  density  was  originally  east  of  the  Earth)  were  to 
"corotate”  with  the  interplanetary  field  lines,  the  fluxes  at  Earth  could 
increase  (3U/3t  >  0)  even  while  3U"/3t  <  0  in  the  corotating  frame  if  G3U/3(j>  < 
3U'/3t  <  0,  since  3U/3t  =  3U'/3t  -  £13U/3$.  We  have  neglected  the  adiabatic 
deceleration  and  radial  gradient  effects  in  (6)  for  the  purpose  of  this 
example.  We  have  documented  many  cases  of  this  "coronal  profile"  effect  (see 
Papers  1,  12,  16  and  17,  for  example),  using  the  field-line  mapping  technique 
described  in  our  previous  Final  Scientific  Report  (Gold  and  Roelof,  AFGL-TR- 
80-100,  March  1980).  If  we  had  anisotropy  measurements,  as  well  as  Intensi¬ 
ties,  the  identification  of  a  rise  as  a  new  injection  is  usually  a  straight¬ 
forward  matter  of  the  presence  of  a  non-equilibrium  anisotropy. 

For  the  decay  phase,  we  cannot  utilize  the  anisotropy  measurements 
without  making  some  assumptions  concerning  the  spatial  gradients  of  the 
intensity,  either  inside  or  outside  of  the  radius  of  the  spacecraft.  The 
simplest  assumption  would  be  that  the  radial  gradient  is  small  in  the  decay 
phase.  Krimigis  and  Roelof  (Study  of  Traveling  Interplanetary  Phenomena/ 1977 , 
Astrophysics  and  Space  Library,  71 ,  D.  Reidel,  1977)  showed  from  the 
University  of  Chicago  simultaneous  measurements  of  11  -  20  MeV  and  30  -  67  MeV 
protons  on  IMP-8  and  Pioneers  10  and  11  that  the  radial  gradients  were  small 
in  the  solar  flare  events  of  November  3,  1973  and  November  5,  1974  between  1 
and  5  AU  when  the  spacecraft  were  in  approximately  the  same  flux  tube.  On  the 
basis  of  conventional  diffusion-convection  theory,  one  expects  a  positive 
radial  gradient  very  late  in  the  event,  but  It  should  also  be  a  small  one, 
with  a  scale  length  >  1  AU  at  such  energies. 


If  we  assume  that  |3(#n  U)/3(0n  r)|  «  1  and  that  the  field  Is 
approximately  spiral  Inside  I  AU,  we  obtain  an  approximate  version  of  (6) 


ir 


1  3* 


Bn  U  +  ft„  C  UJ 


(7) 


where  ftj  and  «2  are  angular  velocities  <  ft,  their  values  being  the  averages  of 
(V  .  )  and  (V  •  V  )  over  the  flux  tube  from  r  to  1  AU,  respectively.  The 
only  quantity  we  can't  measure  directly  with  a  single  spacecraft  is  3  U/3$. 

However,  as  we  have  demonstrated  in  our  earlier  work,  if  the  solar  wind  veloc¬ 
ity  is  varying  sufficiently  to  make  significant  (~  10°)  changes  in  the  connec¬ 
tion  longitude  over  a  few  hours,  it  may  be  possible  to  infer  3U/3 <p  from  the 
event  time  history.  Note  that  3 Bn  U/3t  <  0  and  ft  C  U  >  0  tend  to  cancel  in 
the  decay  phase.  Thus,  if  all  the  terms  on  the  RHS  of  (7)  combine  to  give  the 

observed  anisotropy  £  ,  we  can  assume  that  the  event  is  in  its  decay  phase 

fl  r 

and  that  it  will  continue  to  decay  according  to  the  functional  dependence 
already  exhibited  by  U(t). 


We  tiave  begun  an  examination  of  the  simplifying  assumptions  neces¬ 
sary  to  be  able  to  utilize  the  information  contained  in  S|r»  Multispacecraft 
observations  are  required,  and  this  year  we  finally  received  the  data  which 
enabled  us  to  begin  comparisons  of  IMP  7/8  and  VGR  1/2  particle  measurements. 


To  Illustrate  our  approach,  we  show  in  Figure  1  the  intensity  his¬ 
tories  of  two  channels  responding  to  5  MeV  protons,  P5  on  IMP-8  at  1  AU  and 
CH10  on  VGR1  at  1.5  AU,  the  event  shown  is  the  response  to  the  Ha  importance 
2B  solar  flare  of  0945  -  1108,  22  November  1977  (day  326).  Flare  coordinates 
were  N24 ,  W40  giving  good  connection  to  the  flare  site  at  60°  heliographic 
longitude.  The  flare  was  isolated,  so  the  interplanetary  medium  was  quiet, 
with  V  *  300  km  s~*  so  convection  effects  were  small  until  the  arrival  of  the 
flare  shock.  In  Figure  1,  the  intensities  (which  agreed  quite  well  shortly 
after  launch  of  VGRl)  are  mapped  back  to  their  coronal  connection  longitudes, 
using  hourly  averaged  solar  wind  speeds  from  each  spacecraft.  The  large 
symbols  are  plotted  at  the  beginning  of  each  day,  triangles  for  IMP-8  and 
pluses  for  VGRl . 
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Intensities  are  connected  by  straight  lines  at  the  beginning  of  each 
day.  On  days  327  and  328,  IMP-8  is  higher,  probably  because  of  the  negative 
radial  gradient  expected  near  event  maximum.  However,  by  day  330,  a  rather 
steady  azimuthal  gradient  has  been  established.  There  does  not  seem  to  be  a 
measurable  radial  gradient  after  this  time,  because  the  azimuthal  gradient  is 
not  much  affected  by  significant  changes  in  the  differences  between  the 
convection  longitudes.  The  azimuthal  dependence  can  be  estimated  from  the 
straight  connecting  lines  as  3  0nU/34>  =+  (350)"1,  while  the  residual  temporal 
decay  may  be  deduced  from  the  slow  lowering  of  the  lines,  yielding  3 Bn  U/3t  *- 
(3d)"1.  If  we  substitute  these  values  in  (7)  and  use  V  =  400  km  s  1  for  the 
post-shock  velocity  during  this  decay  phase,  we  find 

£  =  "  ((-0.020)  +  (0.008)  +  (0.005)1  =  +  0.007  (8) 


These  numbers  are  only  preliminary  estimates,  but  they  are  sufficiently 
accurate  to  show  some  consistency  with  (7),  if  /;  =  0.7%. 

We  wanted  to  be  able  to  examine  the  anisotropy  in  some  detail,  so  we 
developed  the  graphical  display  shown  in  Figure  2.  The  eight  gray-shaded 
columns  represent  the  CH10  fluxes  measured  in  the  eight  sectors  scanned  by  the 
LECP  telescope  on  VGR1 .  Large  tic  marks  are  at  one-day  intervals,  and  the 
small  ones  every  6  hours.  The  streaming  directions  run  from  approximately  0° 
on  the  left  through  180°  to  360°  on  the  right.  Protons  flowing  outward  along 
the  nominal  spiral  field  would  appear  at  ~  135°  on  the  plot  (Sector  7).  The 
arrow  points  to  the  direction  for  particles  streaming  away  from  the  sun.  The 
maximum  and  minimum  fluxes  are  indicated  to  the  right  of  the  gray-shade, 

O  __  1 

plotted  as  the  logarithm  in  units  of  (cm  s  sr  MeV)  .  The  gray-shade  assigns 
16  shades  linearly  between  the  minimum  and  maximum  flux  in  each  sector  for 
each  hour.  The  small  squares  each  hour  give  the  direction  of  the  interplane¬ 
tary  field  from  the  MAG  experiment,  projected  upon  the  LECP  scan  plane  (which 
is  tilted  at  this  time  about  10°  from  the  ecliptic).  The  fluxes  for  each 
sector  (and  therefore  the  anisotropies)  have  been  transformed  into  the  local 
E  x  B  frame  for  each  hour. 
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The  strong  Initial  field-aligned  anisotropy  is  evident  during  the 
rise  to  maximum  on  day  326,  and  decreases  thereafter  to  a  maxi mum  ratio  <1.1, 
implying  an  anisotropy  =  0.5%.  Note  how  even  the  very  weak  anisotropies  dur¬ 
ing  the  decay  phase  tend  to  follow  the  outward  field  direction,  so  that  £  > 

0.  Thus,  to  the  accuracy  of  this  preliminary  analysis,  our  use  of  (7),  with 

Its  attendant  assumptions,  appears  to  be  valid  for  this  event.  Assuming  we 
could  obtain  a  reasonable  estimate  of  3  2m  U/3$  from  a  single  spacecraft,  we 
therefore  could  have  used  our  measurement  of  the  parallel  anisotropy  (in  the 
x  B  frame)  to  conclude  that  the  decay  was  primarily  due  to  the  temporal 

term  3  &iU/3t  (the  other  two  being  positive).  Thus  we  would  predict  that  the 

decay  would  proceed  smoothLy,  until  interrupted  either  by  a  change  in  inter¬ 
planetary  conditions,  or  by  a  fresh  injection  of  new  protons. 

Conclusions  and  Recommendations 


The  example  given  above  shows  the  conditions  under  which  we  can  make 
a  straightforward  prediction  for  the  decay  phase  of  an  event.  The  interplane¬ 
tary  medium  is  seldom  in  such  a  simple  state.  For  example,  the  period  one 
month  earlier  of  September  1977  was  strongly  disturbed  by  a  succession  of 
flares  from  the  transit  of  a  major  center  of  activity. 

We  obviously  need  further  diagnostics  of  the  interplanetary  medium 
when  its  configuration  is  more  complex.  We  have  not  yet  analyzed  the  poten¬ 
tial  of  relativistic  electrons  and  lower  energy  (T  <  1  MeV)  ions  in  deducing 

•mch  configurations,  and  their  effect  on  >  l  MeV  proton  fluxes  at  1  AU.  The 

fast  electrons  probe  the  outer  heliosphere,  while  the  slower  <  1  MeV  protons, 

•■/ith  their  stronger  anisotropies,  are  sensitive  indicators  of  the  more 
Immediate  vicinity  (within  ~  0.5  AU)  of  the  spacecraft.  We  believe  the  IMP 
7/8  and  VGR  1/2  data  sets  are  now  ready  for  such  an  investigation,  and  that 
results  to  date  would  justify  a  further  effort  in  this  direction. 
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